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ABSTRACT 

We study mass transfer by Roche lobe overflow in close-in exoplanetary systems. The planet's atmo- 
spheric gas passes through the inner Lagrangian point and flows along a narrow stream, accelerating 
to 100-200 km s~ velocity before forming an accretion disk. We show that the cylinder-shaped accre- 
tion stream can have an area (projected in the plane of the sky) comparable to that of the planet and 
a significant optical depth to spectral line absorption. Such a "transiting cylinder" may produce an 
earlier ingress of the planet transit, as suggested by recent HST observations of the WASP-12 system. 
The asymmetric disk produced by the accretion stream may also lead to time-dependent obscuration 
of the starlight and apparent earlier ingress. We also consider the interaction of the stellar wind with 
the planetary magnetosphere. Since the wind speed is subsonic/sub-Alfvenic and comparable to the 
orbital velocity of the planet, the head of the magnetopause lies eastward relative to the substellar 
line (the line joining the planet and the star). The gas around the magnetopause may, if sufficiently 
compressed, give rise to asymmetric ingress/egress during the planet transit, although more works are 
needed to evaluate this possibility. 

Subject headings: hydrodynamics - planetary systems - stars: individual (WASP-12) - stars: winds, 
outfiows 



1. INTRODUCTION 

The close-in exoplanets (with period less than a few 
days) discovered in radial velocity and transit surveys are 
of great interest, not only because they constrain theo- 
ries of planet formation and evolution, but also because 
they provide a probe of various physical processes that 
are otherwise unimportant in "normal" planets. WASP- 
12b is a transiting exoplanet orbiting extremely close to a 
late-F/early-G star (Af* = I.SSMq, R^ = 1.57i?0, T^s = 
6300 K), with the orbital period P = 1.09 days and or- 
bital semi-major axis a — 0.023 AU = 4.94i?0 = 3.15i?*. 
The planet mass Mp = lAlMj and radius Rp — 1.79Rj, 
as determined by transit observation in optical contin- 
uum (Hebb et al. 2009; Campo et al. 2010). The planet 
is one of the most irradiated exoplanets (with equilibrium 
temperature Tcq — 2500-3000 K) and is highly inflated. 
The fact that the Roche radius (Hill sphere radius) of the 
planet, Rl = a{Mp/3M^y/^ = IMRp, is only slightly 
larger than the Rp derived from optical transit measure- 
ments, suggests that mass loss from the planet is likely 
(Li et al. 2010). The small orbital separation also sug- 
gests that stellar wind may influence the atmosphere the 
planet. 

Recently, Fossati et al. (2010) obtained near-UV trans- 
mission spectroscopy of WASP- 12b with the Cosmic Ori- 
gins Spectrograph on the Hubble Space Telescope. The 
data revealed enhanced transit depths (by about a factor 
of 2) in two wavelength bands, NUVA (2539-2580 1) and 
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NUVC (2770-2811 A), which were attributed to flux at- 
tenuation by absorption lines of metals in the vicinity of 
the planet. Most interestingly, the NUVA data exhibits 
an earlier ingress compared to the transit in optical J,B 
and Z bands, while the egress of the transit occurs at 
about the same time as the optical transit. 

The asymmetric behavior of the ingress/egress in the 
NUVA band relative to the continuum is difficult to un- 
derstand if the absorbing gas surrounding the planet 
arises entirely from an irradiation-driven wind - such a 
wind has been studied extensively in the context of hot 
Jupiter HD 209458b (e.g., Yelle 2004; Tian et al. 2005; 
Garcia Munoz 2007; Murray-Clay et al. 2009): The wind 
is most strongly generated at the planet's dayside, and 
would be distributed on both sides (terminations) of the 
substellar line (the line joining the planet and the star). 
If anything, it would be preferentially on the west side 
because of the planetary rotation, which is almost likely 
synchronized with the orbit (see Schneiter et al. 2007 for 
a simulation). 

In this paper, we consider two possible explanations 
for the asymmetric excess absorption during WASP- 12b 
transit. First, we study the Roche lobe overflow from 
the planet to the parent star (see Li et al. 2010) and the 
associated accretion stream. The stream is asymmetric 
with respect to the line joining the planet and the star. 
We show that the stream has a sky-projected area com- 
parable to the projected planet area, and a sufficiently 
large column density to cause blockage of the star light 
prior to optical ingress. Second, we qualitatively discuss 
the magnetopause produced by the interaction between 
the stellar wind and the planet's magnetosphere. Be- 
cause the planet's orbital velocity, Worb — (GAf^/a)^/^ = 
228 km s~ , is comparable to the stellar wind velocity 
(uu, ~ 100 km s~^ at a = 0.023 AU), the head of the mag- 
netopause lies eastward relative to the substellar point. 
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Again, ingress/egress asymmetry m.ay be produced if the 
gas density in around the magnetopause is sufficiently 
liigfi. 

This paper is organized as fohows. In Sect. 2 we dis- 
cuss Hne absorption by a moving medium with a velocity 
gradient and derive the observational constraints on the 
absorbing gas in WASP-12b. In Sect. 3 we study the 
property of the accretion stream and the obscuration of 
the star light by the stream. Section 4 examines the 
possibility that the gas around the magnetopause may 
absorb the star light. We conclude in Sect. 5. 

2. OBSERVATIONAL CONSTRAINT ON THE EXCESS 
ABSORBING GAS 

We first consider line absorption of stellar radiation by 
the gas around the planet and beyond (e.g., in the accre- 
tion stream). The absorbing gas has a finite temperature 
T and line-of-sight bulk velocity V\\ , both of which (in 
general) depend on the spatial position r. We denote the 
velocity distribution of the column density of the gas (of 
a given species i) by 
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where the integration is along the line of sight, and 
dni/dV\\ is the velocity distribution of the gas density. 

The stellar radiation intensity ft is attenuated to ft = 

ft exp(— r^) after passing through the absorbing gas, 
with the optical depth given by 
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The line absorption cross section can be written as (e.g., 
Rybicki & Lightman 1979) 

<y.{V^^) = .f—<^W^Ml + V^^/c)], (3) 

where / is the oscillator strength of a specific line transi- 
tion, and $ is the Voigt profile centered at 1/0(1 + V||/c) 
(with 1^0 the intrinsic line frequency). When thermal 
broadening dominates, we have $ = l/{^/^TAI/T) at the 
line center, where Avt is the thermal width. Since Avt 
is much less than the bandwidth of transiting observa- 
tions, the Voigt profile can be approximated by a delta 
function. Thus we have 
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where {ne^/meCVo) = (10"i6/4)(Ao/2800l) cm^. 

For WASP-12b, Fossati et al. (2010) found that the 
transit depths in the wavelength regions NUVA (2539- 
2580 1) and NUVC (2770-2811 1) are larger than the 
continuum depth by a factor of two at the 2.5cr level, 
implying the effective radii of order 2.69i?j and 2.66-Rj, 
respectively. The flux depletions in the 41 A window 
of NUVA and NUVC are produced by many absorption 
lines. If we require a single line (of a given atomic species 
i) to significantly deplete the flux (i.e., r^ ^ 1) in the 
wavelength range A A around its intrinsic wavelength Aq, 



then the column density of the gas with velocity around 
V^l = cAA/Ao must satisfy 
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For example, for the Mg II 2800 A lines, / = 0.93, 
and if we take AA = 1 A, we find (diVi/dln V||) ^ 
1.4 X 10^3 cm^^ at V\\ ~ 108 km s^^. For a characteris- 
tic length scale of 3Rp — BARj (see below), this corre- 
sponds to the gas density dni/d\nV\\ ^ 400 cm^"^. Note 
that, since the thermal velocity of the Mg atom is vt — 
(fcr/24mH)i/2 = 1.9(T/104K)1/2 km s'^ such a large 
AA/Ao (= 1/2800, corresponding to V\\ = cAA/Aq = 

108 km s^ ) cannot be due to thermal broadening, but 
must arise from significant bulk motion of the absorbing 
gas along the line of sight. If we choose a smaller AA, 
(i.e., a single line covers a smaller range of wavelength in 
the absorption), the required column density is smaller, 
but then we would require a denser spectrum of absorp- 
tion lines so that the whole 41 A wavelength regions of 
NUVA and NUVC are significantly absorbed. 

Now consider the ingress/egress asymmetry in the 
transit of the NUVA band. Figure 2 of Fossati et 
al. (2010) shows that in optical continuum, the ingress 
starts at the orbital phase (p ~ 0.945, and the center of 
the transit is at cj) — 1. The phase difference of 0.055 
corresponds to a transverse distance (projected on the 
plane of the sky) of 0.925i?^ + Rp ^ IMR^ (for the or- 
bital inclination angle i = 83.1°). In the NUVA band, 
the ingress starts at (^ ~ 0.925 or earlier. The phase 
difference of 0.075 then corresponds to a transverse dis- 
tance of (0.075/0.055)(1.04i?/) = 1.42i?^. This implies 
that the excess NUVA absorber during the ingress must 
extend a distance of at least 0.38i?i 
planetary surface. 

The above constraints for the column density and 
size/location of the line absorber will be used to evaluate 
two possible models for the ingress/egress asymmetry in 
the following sections.. 

3. ROCHE LOBE OVERFLOW AND ACCRETION STREAM 

3.1. Accretion Stream and Disk Formation 

The gas dynamics of Roche lobe overflow in 
semidetached binaries was studied by Lubow & Shu 
(1975,1976). Here we adapt their theory to mass transfer 
in close-in exoplanetary systems like WASP- 12b (cf. Li 
et al. 2010). 

Consider the atmosphere of the planet in quasi- 
hydrostatic equilibrium extending to its Roche lobe. As 
the gas flows through the LI point, it transitions from 
subsonic to supersonic, with the sonic surface close to 
the LI point. The size of the transonic nozzle at LI is 



3.2Rp from the 
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where Cs is the sound speed and fig = 2tt/ P is orbital 
frequency of the binary. The mass transfer rate is related 
to the gas density at LI by 



M~7rpLC.(Ai?L)' 



(7) 
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The steady-state equation of motion of the gas in the 
rotating frame reads 



— = (u- V)u = -V$ 



2r2 X u- 
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where P = pc^ is the gas pressure, and $ is the Roche 
potentiah 
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Here ji = M^,/{M^ + Mp), ri = r = y/JO + Y^, ra 
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+ Y^, and the coordinates {X,Y) are mea- 



sured from the center of the star (see Fig. 1). The LI 
point is determined from d^/dX — for Y ^ 0, giving 
the distance of LI from the center of the planet: 
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for q — Adp/Mi, <C 1. The gas along a streamline 
conserves its Bernoulli integral (assuming constant Cg), 
B = u^/2 -|- $ -|- Cg In p. As the gas leaves LI and ac- 
celerates toward the star, pressure becomes negligible at 
a distance larger than Ai?L, and then the gas particle 
follows a ballistic trajectory. Due to the combined ef- 
fects of Coriolis force and Roche potential force, the gas 
particle leaves the LI region in a narrow stream (with 
width A_Rl), along the direction given by the angle 9s 
(the angle between streamline near LI and the X-axis) 
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with A^a^ [n/rl + (1 - p)lrl]^^. Thus 9^ = 26.7° for 
q = 10"^ (and 28.4° for q -)■ 0). Figure 1 illustrates (in 
the case of g = 10^"^) the stream trajectory for several 
values of initial gas velocities comparable to the sound 
speed [cs/{aQ.o) = 0.02, or c^ ~ 5 km s"^ for the WASP- 
12 system.] 

The crossings of streamlines lead to shock dissipation 
and disk formation. The stream collides with the disk 
at point D (see Fig. 1), whose position (and thus the 
outer radius of the disk r^) is determined by the "no-slip" 
condition, i.e., the (/)-component of the incident stream 
velocity equals the circular velocity at the same point, 
the latter being given by (in the corotating frame) Wdrc — 

(GM/rY^'^—rQ.o (for q<^l). Figure 2 shows the velocity 
of the accretion stream before it collides with the outer 
edge of the disk for the cases oi q — 0.001 and 0.005. We 
see that, for example, for q — 0.001, the stream impacts 
the disk at \Yd\ = 0.24a and r^ ~ 0.71a, with the terminal 
velocity u = O.SSaOo (and Ux reaches O.SSafio). 

As the stream evolves from LI to impact (point D), its 
width W undergoes relatively small variation, due to the 
combined effects of the finite sound speed (which tends 
to make the stream expand) and the enhanced gravity 
closer to the star (which tends to compress the stream). 
Thus, W ~ ARl, and the cross-sectional area of the 
stream remains of order tt^ARl)^. The mean density of 
the stream therefore varies as 
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Fig. 1. — Stream trajectory and disk formation in the orbital 
plane for the case of planet-star mass ratio q = Mp/Af* = 10^'^ 
(WASP-12b parameter). The planet is located at {X, Y) = (-1, 0) 
(in units of a), and rotates counter-clockwise around the star. 
Three streamlines leaving LI are shown, corresponding to the 
initial velocity (in units of the orbital velocity of the planet, 
^/GMJa = 228 km s"!) of {u^,Uy) = (0,0) (solid line), (0.02,0) 
(dashed line) and (0.02,-0.02) (long-dashed line). In the absence 
of a disk, matter would follow these stream trajectories. Self- 
collisions of the streams cause the formation of a disk. The disk 
outer radius is at r^ = 0.71. The stream strikes the disk at point 
D, beyond which the stream does not exist. The dashed straight 
line labeled ti indicates the line of sight when the star light is 
first blocked by the stream (the ingress of the "stream transit"), 
and the line labeled (2 indicates the ingress of the "normal" planet 
transit. A general line of sight is labeled t, with the angle between 
the stream axis and the line of sight denoted by a. 
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Fig. 2. — Stream trajectory and velocity for mass ratio q = 0.001 
(solid lines) and 0.005 (dashed lines). The stream starts at LI and 
terminates when it collides with the outer edge of the disk (the 
vertical dotted line labeled "D"). The star and planet are located 
at (X, Y) = (0,0) and (-1,0) (in units of a), respectively. The 
plotted quantities are u (stream velocity in the rotating frame, in 
units of the orbital velocity of the planet, Q,oa = 228 km s~ for 
WASP 12-b), Ux (the X-component of the stream velocity) and 
\Y\ (the "K-position of the stream; see Fig. 1). 
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for the stream velocity w <; c^. 

Note that the "no-shp" condition Ukely only gives the 
minimum value of the outer disk radius (see Shu & 
Lubow 1981). In reality, tidal forces from the planet set 
the outer truncation of the disk. Numerical simulations 
of the planet-disk interaction in the Type II regime (in 
which the planet is sufficiently massive to open a gap in 
the gas disk; see Lin & Papaloizou 1993) typically find 
that the gap half width is about 0.2-0. 3a (e.g., Lubow 
& D'Angelo 2006; Armitage 2007). Thus we shall adopt 
the outer disk radius r^ ~ 0.7a in this paper. 

The accretion disk formed from the stream has a 
steady-state surface density Sd(r) ~ M /{?>TTVd)^ where 
Vd is the disk viscosity. Adopting the a- viscosity ansatz, 
Vd = ctdHdCds, where Cds is the disk sound speed and 
Hd = Cds/^dif) is the disk thickness, we find 



Sad \CdsJ V^o 



(13) 



where Sl ~ PlW is the gas column density at the LI 
point. Obviously, for a^ ^ 1, we expect S^ 3> Sl. 

3.2. Application to the WASP-12 System 

To estimate the gas density of the accretion stream in 
the WASP-12 system, we assume that the planet has a 
isothermal atmosphere with temperature T ~ 3000 K 
(corresponding to a sound speed of Cs ~ 5 km s~^)(3 
and surface density p{Rp) ^ 5 x lO^'* g cm^'^ (see Li 
et al. 2010). Relative to the planet center, the Roche 
lobe extends to LI at xl — a(g/3)^/'^ — 1.85 Rp. Along 
the y-direction (in the orbital plane), it extends to y^ — 
{2/3)xL = 1.23i?p. Since the Roche potential along the 
planet's y-axis (relative to the planet center) is simply 
— Gm/y-|-constant, one finds p{yL) ~ 3 x 10^^^ g cm^'^ 
(Li et al. 2010). Since the Roche lobe is an equipotential 
surface, this is also close to the density at LI. With the 
nozzle radius ARl ~ 7.5 x 10** km ~ 0.6i?p and a thermal 
outflow velocity c^, we find that the mass transfer rate 
through il is M ~ 2.7 x lO" g s'^ ~ 3 x lO'^Afp yr^^. 

Note that Li et al. (2010) considered mass transfer 
driven by an energy source within the planet's Roche 
lobe. Such an energy source may lie inside the planet, 
maintaining its inflated radius. The suggested energy 
source was tidal dissipation due to an assumed eccen- 
tricity of the orbit (see also Gu et al. 2003). It has since, 
however, become clear that the orbit of WASP-12b is 
highly circular (e — 0.017±Q-Q}f from new radial- velocity 
monitoring; Husnoo et al. 2010), implying that such an 
eccentricity-driven tidal energy source is negligible. Li 
et al. obtained a much higher mass loss rate (about 25 
times larger than our estimate) as they assumed that the 
matter flows from below the Roche lobe over the entire 
surface area of the Roche lobe. In our discussion below, 
we take a conservative approach and adopt our (smaller) 
value of M given above, as we believe it sets the mini- 
mum value of the mass transfer rate. 

The above discussion neglects the possible effect 
of the planet's thermosphere, where photo-ionization 

^ The equilibrium temperature of the planet ranges from 2500 K 
to 3000 K, depending on the albedo and the efficiency of heat 
redistribution in the planetary atmosphere. 



heats the gas temperature to ~ 10"* K (e.g., Murray- 
Clay et al. 2009; Koskinen et al. 2010). Tak- 
ing the photo-ionization cross section crpi ~ 6 x 
10-i8(eo/13.6 eV)-3 cm^ (~ 2 x lO"!* cm^ for photon 
energy eo — 20 eV), we flnd that the column density 
that UV photons can penetrate is 5 x 10^^ cm^^. For 
PL "^ 3 X 10^^^ g cm^'^, this corresponds to a distance of 
2.8 km, much less that AR^. Thus, only the "skin" of 
the accretion stream will be affected by photo-ionization. 
This also shows that, the planet's thermosphere, if any, 
lies outside the planet's Roche lobe. 

In principle, the Roche lobe, the accretion stream and 
the disk can all block the star light. Here we focus on 
the accretion stream since it is asymmetric relative to 
the line joining the star and the planet. To estimate 
the transiting property of the stream, we model it as 
a cylinder with length Sd and radius W — AR^. The 
result of Sect. 3.1 gives Sd = 0.36a (for Xd — —0.668 and 
Yd = —0.238; see Fig. 2). Thus the area of the stream as 
projected in the plane of the sky is SdW sin a ~ GR^ sin a, 
where a is the angle between the line of sight and the 
symmetry axis of the cylinder. Note that a varies during 
the transit (see Fig. 1). Clearly, if the stream is opaque 
to the radiation, it would lead to significant blockage of 
the star light, comparable to that produced by the planet 
itself. 

To evaluate the column density of the stream, we 
need to know how the density varies along the stream. 
As discussed in Sect. 3.1, in the supersonic region, the 
density decreases as u~^ as the stream velocity u in- 
creases. Figure 2 shows that u increases approximately 
as a linear function of s, the distance from LI along 
the stream. Thus we adopt u{s) ~ Cs/F{s), with 
F{s) = 1 for < s < VF and F{s) = W/s for 
W < s < Sd- This gives p ~ plF{s). The hydrogen col- 
umn density in the stream along the line of sight is then 
Nnis) ~ {p/mH){W/ sina) -- lO'^'^ F (s) / sin a cm'^. 
With 1^1 — Cs{s/W) cos a for s > W^, we find 
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for V|| <; 5 cos a km s (where we have used Cs — 
5 km s~^). Thus, as long as the metal abundance satisfies 

Ni/Nh k, 10-8/-nanQ;(V||/100kms"^)^ the stream 
column density would be larger than that required by 
Eq. ([5]). For example, the solar abundance of Mg rela- 
tive to H is 3 X 10-5, and the Mg II 2800 A lines (with 
/ = 0.93) may then produce significant absorption of the 
star light. For Scandium, the solar abundance is about 
10^^, thus stronger lines with / <; 10 or a super solar 
abundance would be needed to produce significant line 
absorption (note that WASP-12a is supersolar by about 
0.3; Hebb et al. 2009). We note that if the gas velocity 
(at a given s) is uniform within the cross-section of the 
stream, then the difference in the gas velocity along a 
given line-of-sight is AV|| ~ Cs cos^ a/ sin a. This is too 
small to produce significant reduction in the observed 
flux. However, it is possible that the 3D gas velocity 
structure in the stream is more complicated (e.g., due 
to the turbulence developed at the surface of the stream 
via Kelvin-Helmholtz instability) and a larger velocity 
difference along the line-of-sight could be produced. 
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The colhsion between the accretion stream and the disk 
produces a hot spot, which may lead to time-dependent 
obscuration of the star hght and possibly early ingress 
(C. Baruteau & D. Lin, 2010, private communication). 
As shown above [see Eq. ([T^]. the accretion disk formed 
from the stream has a significant surface density, and 
can cause obscuration of the star light. Obviously, a per- 
fectly symmetric disk would only give time-independent 
obscuration. However, disk asymmetry (e.g., produced 
by the hot spot) may cause an apparent earlier ingress of 
the planet transit. Numerical simulations will be needed 
to assess these possibilities. 

Our discussion so far has focused on the accretion 
stream due to mass transfer through the inner La- 
grangian point LI. In general, an outward-spiraling 
stream may also form due to mass loss through the 
L2 Lagrangian point (S. Lubow, private communica- 
tion) 0. The L2 point is located at a distance xl2 — 

a {q/3) 1 + (1/3) {q/3) from the center of the 

planet, farther away from the host star than the planet. 
The (Roche) potential difference between LI and L2 is 
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If we assume that the Roche lobe overflow through 
LI does not affect the density of the outer L2 Roche 
lobe, we can estimate the density at L2 as pL2 — 
PLiexp(-2GA'fp/3ac2). For the WASP-12b parameters, 

we find ph2/phi — 0.2 for Cs — 5 km s~^ (see also Gu et 
al. 2003). In reality, pl2/pli will be somewhat smaller 
than the above estimate, and a lower night-side temper- 
ature of the planet will also reduce /9l2- We conclude 
that the spiral stream from L2 is less important than 
the accretion stream. Numerical simulations of the two 
streams would be useful to accurately evaluate their rel- 
ative effects on the planet transit. 

4. MAGNETOPAUSE 

The star WASP-12a resembles the Sun in many as- 
pects, and we will parametrize its wind property using 
the observed solar wind parameters. The fiducial solar 
wind mass flux is M^, = 3 x lO""M0yr-^ At 1 AU, 
the typical solar wind speed is of order 450 km s^ (but 
can change on various timescales from 200 km s^ to 
1000 km s~ ), and the electron density, temperature and 
magnetic field in wind are ~ 7 cm~'^, 1.4 x 10^ K, and 
7 X 10~^ G, respectively (e.g., Hundhausen 1995). Thus, 
the solar wind at 1 AU has a plasma f3 (the thermal pres- 
sure divided by the magnetic pressure) of order unity, and 
the fast magnetosonic Mach number is about 6 (Russell 
& Walker 1995). 

To obtain the stellar wind property at the distance ap- 
propriate to WASP-12b (a = 0.023 AU) requires detailed 
wind models, which we will not attempt here. But it is 
very likely that at such a small distance, the stellar wind 
is subsonic and sub-Alfvenic, and the fast magnetosonic 
Mach number is less than unity (e.g., Mestel & Spruit 
1987; Cranmer 1998). For example, using the simplest 
(isothermal) Parker wind solution (for Mi, = 1.35Mq), 

'' G. Trammell, P. Arras and Z. Li (2010, to be submitted) have 
also recently considered this issue. 



we find that if the wind speed at 1 AU is 500 km s , 
then the sonic point is at 0.032 AU, and the wind ve- 
locity and Mach number at 0.023 AU are 90 km s~ 
and 0.66, respectively. (For a wind speed of 600 km s~^ 
at 1 AU, the corresponding numbers become 0.024 AU, 
148 km s~ and 0.94.) The solar wind temperature ac- 
tually increases with decreasing distance, and this will 
make the Mach number at 0.023 AU even smaller. The 
plasma /? in the wind may reach below 0.1 at such a dis- 
tance, with the wind magnetic field dominated by the 
radial component. 

Given that at 0.023 AU, the stellar wind has a fast 
magnetosonic Mach number less than unity, we do not 
expect bow shock to form when the wind interacts with 
the planet's magnetic field. Instead, near the planet, 
the wind will be "gradually" stopped at the magne- 
topause. Importantly, since the orbital velocity of the 
planet, i>orb = 228 km s~ , is comparable to the wind 
velocity at 0.023 AU, v^ '-^ 100 km s~^, the head of 
the magnetopause is shifted eastward with respect to the 
substellar point by an angle tan~^(worb/i'to) (= 66° for 
v^ = 100 km s"^). 

In general, the location of the magnetopause is de- 
termined by balancing the momentum fiux (P -I- pv'^ + 
B'u/Stt) across the contact surface, where P is the gas 
pressure, v± is the gas velocity perpendicular to the sur- 
face, and i3|| is the magnetic field parallel to the surface 
(e.g., Shu 1992). The standoff distance r^ (measured 
from the center of the planet) at the head of the magne- 
topause is obtained by 



P«;(«Ji+forb)+-fu,+ 



87r 



Pr, 



Bl 



Rz 



(16) 



where Pp is the gas pressure of the planetary magne- 
tosphere, and we have assumed that the planet has a 
dipolar magnetic field with the equatorial strength Bp. 
In the absence of a detailed wind model, we parameterize 
Eq. (HSl) by 



JwPw^w — Jp 
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R'D 
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1 



'^orb 






M- 



(17) 



10 (for 



Vyj ~ 100 km s ^ and the fast magnetosonic Mach num- 
ber M ~ 0.5), and fp ^ 1. With Af„, = Ana^p^jVyj, we 
then have 



Rr, 



2.6 /p'/'(/»i«»iM^i)-^/^B^f , (18) 



where f^i = /lo/lO, Wt„i = w„,/(100km s""^), Mj^i = 
Af„/(3 X lO-"M0yr-i) and Bpi = Bp/{1G). We see 
that the main uncertainty in our estimate of Vm is the 
magnetic field strength of the planet. For example, if 
Bp ~ 5 G, the magnetopause may extends to 4.4i?p. 
Note that r^ is similar to the required extension of the 
absorbing gas to produce the asymmetric ingress/egress 
in WASP-12b (see Sect. 2). This motivates us to consider 
the possibility that the gas (from the stellar wind) around 
Tm may absorb the star light. 
The number density of H in the stellar wind (before 



Lai et al. 



any compression) at a — 0.023 AU is given by 



i^*'^- 

^H 



M„ 



Ana^v^mH 



= 7.7 X 10^ Mwiv~l cm" 



(19) 



As the stellar wind gas is deflected at the magnetopause, 
it may undergo some compression (of order unity) and 
move around the field lines at the velocity comparable to 
Wu, ^ 100 km s^^. Whether this gas can provide enough 
line absorption to explain the observation depends on 
the physical states of the gas. Given the low gas density 
in Eq. (J19p . we consider the most optimistic scenario: If 
the gas around r^ is able to cool to 10'' K (the equilib- 
rium temperature between photo-ionization heating and 
coUisional cooling) from 10^ K (the original stellar wind 
temperature) while maintaining pressure equilibrium, it 
may be undergo compression by a factor of 100. Thus 
the characteristic hydrogen column density is 

Nh ^ lOOngV™ ^ SxlO^'f^^if-^vZlMl.y^'Bli'cm- 

(20) 
Comparing this with Eq. ([5]), we see that, under this opti- 
mistic condition, sufficient metal column density may be 
attained when ni/nn ^ 5 x 10^^, which is of the same 
order as the stellar abundance of Magnesium Q. Obvi- 
ously, a larger M^ would result in a higher gas column 
density. 

It thus seems that the compressed gas at the magne- 
topause might just be able to cause enhanced absorp- 
tion before ingress. The interaction of the stellar wind 
with planetary magnetospheres is, however, much more 
complex than discussed above (e.g., Kivelson & Russell 
1995). Our estimate in Eq. ([^0]) suggests that such inter- 
action may produce observable effects. Further study on 
this subject in the context of hot Jupiters is worthwhile 
(cf. Preusse et al. 2007; Ekenback et al. 2010). 

5. DISCUSSION 

Motivated by the recent near-UV spectral observations 
of the WASP-12 planetary system (Fossati et al. 2010), 
which showed an earlier ingress of the planet's transit 



compared to the optical continuum, we have studied two 
possible explanations for the ingress/egress asymmetry. 
The first involves mass transfer through Roche lobe over- 
flow, which results in an elongated accretion stream in 
which gas flows from the LI Lagrangian point toward the 
star. Our analysis of the geometric and velocity structure 
of the stream shows that it may indeed provide asymmet- 
ric obscuration of the star light during the planet tran- 
sit. In addition, it is also possible that an asymmetry in 
the accretion disk, caused by the impact of the accretion 
stream, produces an apparent earlier ingress. Another 
possibility involves the magnetopause separating the stel- 
lar wind and the planetary magnetosphere. Because of 
the planet's orbital motion, the head magnetopause lies 
eastward relative to the substellar point. We suggest that 
line absorption by the gas around the magnetopause may 
also explain the asymmetric ingress/egress behavior, al- 
though more theoretical works are needed to understand 
the property of the absorbing gas around the magne- 
topause. 

The two possibilities studied in this paper may be 
distinguished by the fact that in the accretion stream, 
the gas falls away from the observers, producing red- 
shifted absorption, while the flow around the magneto- 
sphere tends to move toward the observer, thus produc- 
ing blueshifted absorption. 

Finally, we note the observed behaviors reported by 
Fossati et al are only marginally significant (at most 3a 
effects). Thus, more observations will be useful. Nev- 
ertheless, the physical processes discussed in this paper 
are quite general and therefore may be applicable in other 
close-in exoplanetary systems. 
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"^ What is of interest here is the metal abundance of the stellar 
wind, which is unknown. In the case of the Sun, it is known that 



the solar wind metal abundance can be larger than the abundance 
at the solar atmosphere (e.g., Bochsler & Geiss 1989). 



